Stabilizer quantum error correction toolbox for superconducting qubits 
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We present a general protocol for stabilizer measurements in a system of TV superconducting 
qubits. We assume always-on, equal dispersive couplings X to a single mode of a high-Q microwave 
resonator in the ultra-strong dispersive limit denned by X > l/T2,/t, where T2 is the qubit coher- 
ence time and ft is the cavity linewidth. In this limit, we show how to map the two eigenvalues 
of an arbitrary weight M < N Pauli operator, onto two quasi-orthogonal coherent states of the 
cavity. Together with a fast cavity readout (T raea s <C this enables the efficient measurement 

of stabilizer operators. 



Several milestones on the road to quantum comput- 
ing with superconducting circuits have recently been 
reached, such as the experimental demonstrations of 
Shor's factorization algorithm [1] and of rudimentary 
quantum error correction (QEC) [2]. As the resources 
required for more complete QEC protocols get within ex- 
perimental reach, it is desirable to develop a toolbox suf- 
ficiently versatile to allow the implementation of a wide 
class of codes. Most QEC codes can be described con- 
cisely in the stabilizer formalism of Gottcsman [3]. In 
this framework a QEC code is defined by the subspace 
spanned by the eigenstates with eigenvalue +1 of a set of 
commuting multi-qubit Pauli operators called stabilizers. 
Error detection is achieved by measuring the stabilizers 
of the code; the syndrome of an error being a sign flip of 
stabilizers. Correction in turn, can be performed when 
the collection of syndromes contains enough information 
to identify the location and type of the error. The ability 
to measure arbitrary multi-qubit Pauli operators would 
thus allow a direct realization of stabilizer QEC codes. 

Recent dramatic progress in coherence times for su- 
perconducting qubits coupled to 3D cavities [4-6], has 
inspired new theoretical ideas for mapping qubit states 
to Schrodinger cat cavity states [7, 8]. Following, these 
developments, we show in this work, how to realize stabi- 
lizer measurements with superconducting qubits coupled 
off-resonantly to a common mode of a high-Q cavity. Ad- 
vantages of our approach over previous proposals [9, 10] 
include a simpler hardware and the ability to selectively 
address an arbitrary subset of qubits, without the need 
for tunable couplings. Besides its relevance for stabilizer 
QEC, our approach paves the way for the preparation of 
arbitrary cluster states for measurement-based quantum 
computation [11] and for the digital simulation of quan- 
tum lattice models with non-trivial topology such as the 
toric and surface codes [12, 13], similar to a recent work 
on trapped ions [14]. 

Central to our proposal is the off-resonant coupling 
between a superconducting qubit and a single mode of 
a microwave resonator [15] described by the dispersive 
Hamiltonian i/disp = Xcr z a { a, where cr z = \e)(e\ — \g}(g\ 
is the Pauli matrix in the computational basis {\g), |e)} 
of the qubit and a (a^) denotes the photon annihilation 



(creation) operator of the cavity mode. This coupling 
describes a qubit-state-dependent frequency shift ±X of 
the cavity, or equivalently a photon-number-dependent 
frequency shift 2nX of the qubit. In the weakly dis- 
persive regime \jTi,n ~ 2X, where k is the bare cav- 



ity linewidth and T 2 



(27i)- 



Ta, is the qubit co- 



herence time composed of relaxation 1/T\ and pure de- 
phasing T^, this interaction has been used for a qubit- 
rcadout by measuring the phase-shift of transmitted or 
reflected microwaves [15]. In the strong dispersive regime 
l/T-2, k < 2X, the same interaction can be used to resolve 
different photon numbers in the qubit spectrum [16]. In 
this work, we are interested in the ultra-strong dispersive 
regime, where k,1/T2 <C X. In this regime, it becomes 
possible, by using microwave pulses narrow in frequency 
compared with X but large compared with k and I/T2, 
to perform qubit rotations conditioned on the number 
of photons in the cavity, or cavity displacements condi- 
tioned on the state of the qubit [7, 8]. Here we show 
how to take advantage of this novel regime to encode the 
two eigenvalues of an arbitrary multi-qubit Pauli opera- 
tor onto quasi-classical oscillations of light that differ in 
phase by ir. 

Although our scheme is applicable to other types of 
superconducting qubits, for clarity we will frame our dis- 
cussion around the specific case of transmon qubits. A 
transmon qubit [4, 17] is formed by a superconducting 
dipole-antenna with a Josephson junction at its center 
with Josephson energy Ej ^> Ec = e 2 /(2Cs), where 
Cs represents the total capacitance between the an- 
tenna pads. Neglecting charge-dispersion effects, which 
are suppressed exponentially in Ej/Ec [17], the low- 
energy spectrum of an isolated transmon is well approx- 
imated by that of an anharmonic oscillator with fre- 
quency LJ01 ~ \/oEjEc — Ec and weak anharmonicity 
uqi — U12 ~ Ec *C u>oi- In state-of-the-art realizations, 
the qubit linewidth I/T2 is close to being limited by re- 
laxation [4, 18-20]. In this work we are interested in 
a setup such as depicted in Fig. 1 (a), where TV trans- 
mons are coupled dispersively to a microwave field inside 
a 3D cavity. For control and readout purposes, an ancilla 
qubit, is further dispersively coupled to both the high-Q 
cavity containing the N qubits with Xa 3> AX and to a 
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FIG. 1. (Color online) (a) 2D array of N transmon qubits in 
a 3D cavity. The ancilla qubit and the upper cavity are used 
for readout /reset and manipulation purposes. Cavity (b) and 
qubit (c) spectra in the ultra-strong dispersive regime. 

low-Q (readout) cavity, similar to the setup used in [6]. 
We assume that both the ancilla qubit and the readout 
cavity remain in the ground state, except during readout 
and manipulation. Thus omitting, for now, the corre- 
sponding degrees of freedom, we model this system in an 
appropriately rotating frame (see supplemental material 
for details), by the effective Hamiltonian 

N 

H = tr I a) a - Ka^a^aa . (1) 

i=l 

The transmons are treated here as two-level systems as- 
suming their anharmonicity is larger than their lincwidth 
(i.e. Ec > I/I2). Furthermore, we neglect the cavity- 
mediated qubit-qubit interaction, assuming the qubits to 
be sufficiently detuned from each other. Crucially, we 
require the N qubits inside the cavity to have equal dis- 
persive coupling X with the cavity mode. In practice this 
may require some fine-tuning of the individual qubit tran- 
sition frequencies and coupling strengths. The frequen- 
cies could be tuned for example by using split- Josephson 
junctions for each transmon and flux-bias lines to control 
their effective Josephson energy [21], while the coupling 
strengths could be adjusted through the geometry and 
placement of the transmons in the cavity. The second 
term on the rhs of Eq. (17) accounts for the (negative) 
qubit-induced anharmonicity of the cavity [22, 23]. In 
the weak and strong dispersive regimes, this term can 
usually be neglected as K -C X. We find that in the 
ultra-strong dispersive regime it is necessary to account 
for its leading order effect. We next show how to encode 
the parity Z$ N = ® i=1 of of an iV-qubit state \ip) N 
onto two quasi-orthogonal coherent states of the cavity 
differing in phase by w. 



Parity encoding. Suppose the system is initially pre- 
pared in the product state \^) t=0 — \a)\ip) N , where \a) is 
a coherent state of the cavity with amplitude a. Making 
use of the identity exp[— i(n/2) Y^,iLi °f] = ( — ^) N ^Sni 
one can show that under the action of (17), at time 
T = n/(2X), the state becomes 

|*) r = U K (\a N )P+ N + \-a N )Pg N ) \yj) N , (2) 

where P s = (1 ± Zs N )/2 are the projectors onto the 
even (+) and odd (— ) qubit parity subspaces as measured 
by the ±1 eigenvalues of the multi-qubit Pauli operator 
Zs and a N = (—i) N a. Note that the self-Kerr term is 
qubit-independent [24] and conserves the photon number 
a^a. Because it commutes with the dispersive term, its 
effect factors out and is captured in (19) by the unitary 
operator Uk = exp[iirK/(2X)a'a'aa]. For weak nonlin- 
earity such that nK -C X, the leading order effect of Uk 
acting on the coherent states \±a N ) is a rotation of the 
mean amplitude by an angle A0 n i = ttuK/X with the 
mean photon number n = \a\ 2 . To leading order in K/X, 
the state (19) is thus well approximated by 

m T = \a N )P+ N \i>) N + \-a N )Ps N W) N , (3) 

with c\n = a N e l ^ nl . The sub-leading order ef- 
fect is a damping of the mean amplitude by a factor 
cxp(— A^[/(2n)) [25]. We emphasize that in the ultra- 
strong dispersive regime k/X < 1 considered here, pho- 
ton decay only weakly damps the amplitude of the co- 
herent states in (12) by a factor exp(— kT/2) ps 1 — 
(7t/4)(k/X). Ignoring these small effects, we thus see that 
the dispersive interaction can be used to encode the par- 
ity of the multi-qubit state onto two coherent states of 
the cavity differing in phase by 71". 

Subset selectivity. Often stabilizers are defined on sub- 
sets of qubits. Selectivity to M < N qubits, labeled by 
the set S M C Sn = {1, . . . , N}, can be achieved as fol- 
lows. Consider the identity 

where Us(t) = exp(—itXa^ aY^ies a X) anc ^ = 
Sjy\S M is the complement set to S M . Eq. (24) can be 
easily shown using a x a z a x = —cr z . Thus, by splitting 
the dispersive evolution of all N qubits into two equal 
halves and interspersing it with bit-flip operations on the 
qubits in S\ il we can effectively echo away the contribu- 
tion of the latter to the total "magnetization" and im- 
plement the dispersive evolution Us (t) of the qubits 
in S M alone. Acting on an initial state of the form 
\a)\ip) N , Us M (T = 7r/(2X)) then encodes the subset- 
parity Zs M = ®ies M °f onto tnc state of the cavity 
as explained above. 

Physically, the initial unconditional cavity displace- 
ment and bit-flips can be implemented via fast microwave 
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pulses (see supplemental material). Because of the dis- 
persive interaction, the qubit transition frequency of 
the z-th qubit splits into a ladder of frequencies w™ = 
lo^l + 2nX, corresponding to different photon numbers in 
the cavity (Fig. 1 (c)). The latter are Poisson distributed 
and peaked around n. Hence, to best approximate an un- 
conditional rotation of the i-th qubit, the pulse must be 
centered at the frequency LJg i +2nX and have a frequency- 
width large compared with 2\ffiX. For n > 1/tt the du- 
ration of such a 7r-pulsc is thus T x ^ l/(2^/Ex) < T. 
Similarly, the initial coherent state of the cavity \a) can 
be prepared from the vacuum by driving the cavity at the 
frequency uj c — Xa with a pulse of area a and a frequency- 
width large compared with 2NX; the maximal frequency 
spread of a cavity dispersively coupled with strength X 
to N qubits (Fig. 1 (b)). Again the duration T<i of this 
pulse is short since T d < 1/(2NX) < T. The total dura- 
tion of the encoding is thus dominated by the dispersive 
evolution time T — it/ (2X). 




FIG. 2. (Color online) (a) Quantum circuit diagram for en- 
coding the parity of qubits 2 and 4 (full (blue) lines). D a 
represents the displacement operation, ir x a single-qubit n- 
pulse and Txn a free dispersive evolution of duration 7r/(4x). 
(b) Numerical simulation of the evolution of the Q-function 
of the cavity [25], with an initial qubit state = (\ggge) + 
\gg e g) + \ eee g))/VS. Dissipation from photon loss at a rate 
k/(27t) = 10 KHz and qubit decoherence with Ti = T2 = 20 /is 
are included as well as a finite displacement and 7r-pulse du- 
ration of Ins. Other parameters are: a = 2, x/(27r) = 5 MHz 
and K/(2n) — 80 KHz. The self-Kerr term leads to an addi- 
tional phase rotation A<j> n \ — 2KnAt, where At = 50.3 ns is 
the total duration of the encoding. Taking this rotation into 
account, we obtain a fidelity of F = 98% to the ideal target 

State l*)idcal = |a2> J P { + 2 , 4 }IV'>4 + h«2>P { 2,4 } IV>>4' 

Fig. 2 shows the quantum circuit diagram for the 
subset-parity encoding of M = 2 qubits out of N = 4 
qubits along with the results of a numerical simula- 
tion with square cavity and qubit pulses of duration 
T v = T d = Ins, dispersive shift X/(2n) = 5 MHz and 
coherent state amplitude a = 2. The initial four-qubit 
state is \tp) i = (\ggge) + \gg e g) + |eee<j))/v3 and the 
simulation includes both dissipation from photon loss at 
a rate k/(2tt) = 10 KHz and qubit decoherence with 



Tx = T<x = 20 /US (see supplemental material). Fur- 
thermore we take into account a cavity nonlinearity of 
K/(2ir) = 80 KHz. Accounting for the phase rotation 
induced by the latter, and shortening the dispersive evo- 
lution times to compensate for the finite duration of the 
7r-pulses, we find a fidelity (overlap with the ideal tar- 
get state |tf) idcal = |a 2 )P+ i} \ip) 4 + \-a 2 ) P^ A} \ip) 4 ) of 
98%. So far we have shown how to encode the parity 
of a weight M < N cr z Pauli operator Z3 . By apply- 
ing single-qubit rotations to individual qubits before and 
after the encoding one may in fact encode the parity of 
an arbitrary weight M Pauli operator Qs M = &)ies T *' 
with Ti G {erf , er^ , of}. Next we explain how the en- 
coded parity can be measured. 

Parity readout. The encoded state is of the form 

I*)t = \ & m) P S m \^)n + \-^m) P s u \^)n^ where in gen- 
eral P| M = (1 ± Qs M )/2 and a M = (-i) M e iA ^a. The 
overlap between the two cavity states, (o M \ — a M ) = 
cxp(— 2|a| 2 ), is independent of K and M. For sufficiently 
large amplitude \a\ of the initial displacement, these two 
states are distinguishable in principle and a measurement 
of the cavity state is equivalent to a multi-qubit parity 
measurement. A fast readout of the cavity state with 
T m cas <C 1/k, may be achieved by lowering the Q-factor 
of the cavity containing the qubits (k — > k' k), as 
recently demonstrated [26]. This Q-switching adversely 
affects the lifetime of the qubits via the Purccll effect. 
However, the latter is expected to be weak as long as 
k' <C X. Alternatively, the cavity state can be mapped 
onto the ancilla qubit, which can subsequently be mea- 
sured through standard homodyne measurement of the 
low-Q readout cavity. The mapping is achieved physi- 
cally in three steps. First, the high-Q cavity field is dis- 
placed unconditionally by a M . To a good approximation, 
this maps the encoded state onto 

D 5m |*> t = \2a M )P+ u \^) N + \0)P Sm H') n ■ (5) 

The second step consists in performing a 7r-pulse on the 
ancilla qubit, which so far was in its ground state, con- 
ditioned on the cavity being in the vacuum state. As 
first proposed in [7] and demonstrated in [6], this can 
be achieved by applying a pulse centered on the bare 
ancilla qubit transition frequency, which is narrow in fre- 
quency compared with 8nX A (the additional factor of 
4 is due to the twice as large amplitude of the cavity 
state in the first term on the rhs of Eq. (28)). Because 
X A 3> NX, a pulse duration Ta can be chosen such that 
1/(2Xa) < T A < 1/(2JVX). For n > 1/4 the first in- 
equality guarantees conditionality while the second one 
allows us to neglect the dispersive evolution during this 
operation. The state then becomes approximately 

\2a M ) (P+ A/ \*P) N ) \g) A + |0) (P S " M W) N ) \e) A . (6) 

In the third and final step, a displacement of —2a M is 
performed on the cavity conditioned on the ancilla qubit 
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being in the ground state. This is achieved with a pulse 
centered at frequency u c — Xa, with a frequency- width 
small compared with 2X A but large compared with 2NX 
(Fig. 1 (b)). Neglecting again the dispersive evolution 
during this step, the state finally becomes 



io) {p£ M m N Ma 



'N 



(7) 



Note that the state (25) is now stationary with respect to 
the dispersive interaction, there being no photons in the 
cavity. Reading out the state of the ancilla qubit amounts 
to measuring Qs M - After the measurement, the ancilla 
qubit may be reset to the ground state efficiently via the 
readout cavity [27]. 

Simulated time evolution. Apart from enabling a fast 
parity measurement, the ancilla qubit can also be used 
to simulate the "time evolution" exp[— i8Qs M ] under the 
action of an arbitrary Pauli operator Qs M ■ This is useful 
e.g. for manipulating a logical qubit state encoded in a 
stabilizer subspace in which case Qs M i s taken to be a 
logical qubit Pauli operator. Assuming that arbitrary 
rotations can be performed on the ancilla qubit in a time 
short compared with 1/(2NX), we start with Eq. (28) and 
perform a 27r-rotation of the ancilla around an axis n = 
(s'm(ip), 0, cos(if)) withy) = arccos(l — 28/tt), conditioned 
on the cavity being in the vacuum state. As a result, the 
odd-parity component of (28) acquires a phase 29 (half 
the solid angle subtended by the ancilla state) and the 
state becomes 



\2& M )P+ M \iP) N + e™\0)Ps M \^) N . 



(8) 



Note that in (31) we omitted the state of the ancilla 
qubit, since it starts and ends in the ground state and 
thus factors out. Finally the encoding is undone by ap- 
plying the encoding protocol with a replaced by — a M . 
This disentangles the cavity and qubits and after uncon- 
ditionally displacing the cavity back to the vacuum, tak- 
ing into account an additional nonlinear phase acquired 
during the decoding, the W-qubit state becomes 

PU^N + e 2ie PsJi')N = e W e- ieQs "\i') N , (9) 

which up to an unimportant global phase factor, repre- 
sents the action of the desired unitary. 

Application. To test the feasibility of the above 
protocols, we simulated the preparation of a log- 
ical qubit state of the four-qubit erasure channel 
code [28]. The stabilizer generators of this code are 
S = {Z1Z2, Z3Z4, X1X2X3X4}, where we have switched 
to the standard notation Xi and Zi for the Pauli oper- 
ators of qubit i. The code-space is spanned by the two 
+1 cigenstates of the stabilizers: 

\0) = l(\gg) + \ee))(\gg) + \ee)) (10) 

\T) = \{\gg)-\ee)){\gg)-\ee)). (11) 



The logical qubit Pauli operators are Z = X\X 2 = X3X4 
and X = Z X Z 3 = Z 2 Z 4 = Z X Z± = Z 2 Z Z . Because of the 
redundancy of the logical operators, this code protects 
a logical qubit state a|0) + /3|l) from the loss (i.e. ar- 
bitrary error) of a known qubit [28]. Here we prepare 
the logical qubit state = cxp[— z(7r/8)A]|0) as fol- 
lows, (i) We start with the fully polarized four-qubit 
state \gggg), which is already a +1 eigenstate of the Z 
stabilizers, (ii) Using the encoding protocol, we measure 
the logical Z operator X\X 2 - If we obtain —1, we apply 
Z\. (Hi) We measure the stabilizer X1X2X3X4. If we ob- 
tain — 1 we apply Z 3 . In the absence of qubit- frequency 
tunability, Zi and Z 3 can be implemented, by combin- 
ing two Hadamard rotations and a 7r-pulse [15]. (iv) We 
reset the cavity to the vacuum. These four steps pre- 
pare the logical state 1 0} . We next use the ancilla (with 
lOiVX) to implement the rotation as described 



*a 

above with Qs 4 = X and 9 = 



Fig. 3 shows the ob- 
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FIG. 3. (Color online) Fidelity of the prepared state to 
cxp[— i (71-/8) A"] 1 0). The (red) dashed curve represents the 
boundary of the inequality K > x 2 /(4a 9 ) which relates the 
self-Kerr K with the dispersive shift x an d the single-qubit 
anharmonicity a q [22]. Here a,/(27r) = 200 MHz. 

tained fidelity to the ideal target state exp[— z(-7r/8)A] |fj) 
as a function of X and K in units of k/(2tt) — 10 KHz and 
for Ti = T 2 = 20 jis. The total duration T of the state 
preparation is shown on the upper x-axis in units of l/n. 
Focusing on the line K = 0, when X is small T is large and 
dominated by the dispersive evolutions and the fidelity 
is limited by a combination of qubit decoherence, pho- 
ton loss and faulty conditional ancilla rotation. The fi- 
delity then increases with increasing X, which reduces the 
preparation time T and hence the effects of decoherence 
and improves the fidelity of the conditional ancilla rota- 
tion. It reaches a maximum in a regime where both con- 
ditional and unconditional operations can be performed 
with high fidelity. A further increase in X degrades the 
unconditional cavity displacement and qubit rotations (a 
1 ns pulse corresponds to a width of ~ 160 MHz) and the 
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fidelity drops. The cavity nonlincarity K and the dis- 
persive shift X are in fact not independent, but rather 
related via the single-qubit anharmonicity a q , by the in- 
equality K > X 2 /(4a g ) [22], which is shown as a dashed 
(red) curve in Fig. 3 for a g /(27r) = 200 MHz. 

Conclusion. We have proposed a protocol based on 
an echo technique to realize arbitrary stabilizer measure- 
ments with superconducting qubits in the ultra-strong 
dispersive regime. Numerical simulations, including de- 
coherence and non-linear effects, show encouragingly 
high fidelities. 
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In these notes, we derive the effective Hamiltonian of Eq. (1) in the main text. We calculate 
finite pulse-width effects for the unconditional displacement and single-qubit rotations and discuss 
the constraints imposed on the drive and system spectrum. We describe a modification of the parity 
measurement protocol for stabilizer pumping. We briefly discuss our treatment of decoherence and 
present the state tomography of the encoded logical qubit state. Finally we conclude with a discussion 
of possible applications and extensions. 



DERIVATION OF THE EFFECTIVE MODEL OF EQ. (1) 

In terms of its two conjugate variables ip and n, which represent the superconducting phase difference and the 
number of Cooper pairs transferred across the Josephson junction (JJ), a transmon qubit is conventionally described 
by the Hamiltonian [17] 



H = AE c {n - n g ) 2 + Ej{1 - cos( V )) = AE c {n - n g f + ^-<p 2 + KiM • (12) 



In the second equality we have explicitly separated the harmonic and anharmonic parts of the Josephson potential 
Ej(l — cos(y)) = (p 2 Ej/2 + V n \((p). For Ej ^> Eq, the "phase particle" essentially undergoes quasi-harmonic 
oscillations around a local minimum of the cosine potential, with rare 27r quantum phase-slip events, which lead to 
a weak (exponentially suppressed in Ej/Eq) dependence of the energy levels on the offset charge n g [17]. In the 
following we will neglect the latter and set n g = 0. It is convenient to diagonalize the harmonic part by introducing 
the bare mode operator ao = y/ Ej / (2fruo) <P + 2iy Ec /(Hu>o) n and frequency ftwo = y/8EjEc- Then expanding the 
non-linear potential to leading order (<p 4 ) one obtains 

H = hw {a\a + 1/2) - ^(a + aj) 4 + 0(<p 6 ) (13) 



Note that the coefficient of a given term in the expansion is smaller than the previous one by a factor ~ y/Ec/Ej. We 
are primarily interested in low-lying energy levels, for which the wavcfunctions are localized in phase i.e. y/ (f 2 ) "C 7T. 
Hence we consider here only the leading order i^ 4 non-linearity. 

The coupling of a transmon with the electromagnetic field of a resonator is typically described in the dipole 
approximation [30] by a term ~ n • E where E = ^ k iek(bk ~ bt) is the quantized electric field of the (multi-mode) 
resonator written in terms of its eigenmodes b^. In a straightforward generalization of the above, the Hamiltonian of 
7V q b transmons coupled to multiple harmonic modes can be written in the general form (suppressing from now on the 
zero-point energies) 



H = J2 fafalat + ]T hufifr +5> J -K - a+X&j - b]) + J2 V®{ Vi ) with 9ij = | 

i=l j ij i=l 



qb 



(14) 



Here V^\<pi) = Ej[(l — cos(y>j)) — ip 2 /2}. As before, it is convenient to diagonalize the harmonic part. As 
shown in [22] , this diagonalization can be achieved from the knowledge of the iV q b X iV q b impedance matrix Z of the 
corresponding linearized circuit with each one of the iV q b JJs representing a port. Denoting the normal modes by c p 
and their frequencies by ui p , one obtains explicitly 



Here (f>o = h/(2e) is the reduced flux quantum and we have introduced the admittance at port k defined as Yfe = 1/Zj-k- 
The choice of the port k is arbitrary and reflects a particular basis choice. The total number of modes is denoted by 
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M. The mode frequencies uj p are given by the zeros of the imaginary part of the admittance, i.e. ImYk(u> p ) = (note 
that in the dissipationless case, the elements of Z are purely imaginary). 

Keeping the leading order ip 4 non-linearity and normal-ordering the latter, the Hamiltonian of the system then 
takes the form 

H± = ^2 Upclcp - ^2 Ipp' ( 2c I c p' + 4 c p' + c p c p' ) ( 16 ) 
p pp' 

^ ] fipp'qq' {^> c pCpiCqCqi + Ac^C^, C^Cqi + ACpCpi CqCqi + CpCpi CqCqi + CpC^C^C^j , 
pp'qq' 



with 



pp'qq' 



— Rn 



^7cff ^7cff <7cff <7cff \ , 
^kp^kp'^kq^kq' 2-^1 



1 ^j' Zik(uJp) Zik(uJp>) Zik(bJ q ) Zik(iO q i) 



6 Zkk{Up) Zkk{u p i) Zkk{Uq) Zkk(Wq') 



(17) 



where i?o = h/e 2 and 7 PP ' = §^2 q Pqqpp' ■ Eqs. (16) and (17) together with a model for the impedance matrix, allow 
one to design the circuit Hamiltonian [22, 23]. Importantly, the above derivation provides the functional dependence 
among the parameters ui p , "/ p > and (3pp> qq i. By applying a generalized RWA in the dispersive limit of large detuning 
between the modes [31], we may neglect rapidly- rotating and non-diagonal terms and thus consider the simplified 
form 

H cS = ^ UpCpCp + 2 S ^pq c l c \ c i c P ( 18 ) 
p pi 



where we have defined X pq = — 24/3 pggp for q ^ p and X pp = — 12/3 pppp . The frequencies have been renormalizcd 
as oj'p — uip — 2 r y P p. The second term on the rhs of (18) describes self-Kerr (X pp ) and cross-Kerr (X qp with g ^ p) 
interactions between the modes. From the Cauchy-Schwarz inequality and Eq. (17), it can be shown that X 2 p < 
AXqqXpp] equality holding in the singlc-qubit case. In the dispersive regime of interest here, the hybridization of 
qubit and cavity modes is small and one can unambiguously identify the qubit-like modes as those with the strongest 
anharmonicitics. Assuming the latter are much larger than the linewidth and as long as the drive frequencies used 
are different from transitions between the first and second excited states, the qubit Hilbert space can be truncated to 
the lowest two levels. From (18) we then obtain, denoting the cavity modes with a, and suppressing constant terms 



N, 



q b qb 



(19) 



where i-Q = Xu and we have kept only the dispersive and qubit-cavity cross-Kerr shifts. We note that in addition 
to these interaction terms there exist in general also qubit-qubit (~ cftr?, i ^ j) and cavity-cavity (~ a\aia}-aj^ 
i 7^ j) cross-Kerr shifts. However in the regime of interest here, the coefficients of these terms are small and will 
consequently be neglected. 

A qubit-induccd cavity non-linearity has previously been derived within the dispersive approximation of a multi- 
level generalization of the Tavis-Cummings model [32]. In addition to the qubit independent term discussed here, such 
a treatment also yields qubit dependent non-linearities. In our derivation, qubit dependent self-Kerr terms appear 
only at order ip 6 and higher and are paramctrically small for large ratios Ej/Ec- 

We can now apply (19) to our system depicted in Fig. 1 of the main text. From the above derivation, we see that 
in principle all qubits (N active qubits and 1 ancilla qubit) couple to multiple modes of both cavities. However, in 
the weakly coupled dispersive limit for the readout cavity and appropriately chosen parameters, we may neglect the 
coupling of the active qubits with the readout (low-Q) cavity mode and higher modes of the active (high-Q) cavity. 
Similarly the ancilla approximately only couples to a single mode in each cavity. Choosing the dispersive shifts of the 
active qubits to be all equal to X and denoting the dispersive shift of the ancilla with the active cavity mode a by 
Xa NX and with the readout cavity mode b by X' A we obtain 

N qb N 

H a =^2 ~Y CT t + ^1*1 + a-a^a + X A ^a^a - Ka)a)aa + ^readout , (20) 

4=1 1=1 
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with ifroadout = w' c tf b + X' A a^b^b — K'b^b^bb. Finally, moving to a frame rotating with qubits and cavity by the 
unitary transformation U(t) = exp —it ( /2)of + (ui c — ^A)a)a\ yields Eq. (1) of the main text: 



N 



Hq = X cr^a^a — Ka'a'aa . 



(21) 



This model is a valid when the ancilla qubit and readout cavity remain in their ground state, in which case i? rca dout = 0. 



FINITE PULSE- WIDTH EFFECTS ON UNCONDITIONAL OPERATIONS 



Having motivated the form of our Hamiltonian (21) [Eq. (1) of the main text], we next use it to investigate the effects 
due to finite pulse- width when driving the system with realistic pulses. To distinguish these effects from the nonlinear 
effects we here set K = 0. In the following discussion, we concentrate on the unconditional operations necessary for 
the parity encoding protocol and neglect the ancilla qubit and readout cavity. We further assume independent drives 
for the qubits and the cavity mode. For the latter to be a valid approximation, the qubits need to be sufficiently 
detuned from each other and from the cavity. In the rotating frame of Eq. (21), a cavity drive at frequency iVd is then 
described by 

H c d (t) = e (t)(e it ( u '- XA - Ud )a J t + h.c). (22) 

In the same rotating frame, a simultaneous drive of the qubits in a set S £ Sn = {1, • • • ,N} at frequencies tof for 
i £ S, is described by a term 

Hf(t) = £ ^( CT +e^f + h.c). (23) 



Unconditional displacement 

At the start of the encoding, we need to be able to displace the cavity from the vacuum into a coherent state with 
amplitude a independently of the state of the qubits. We choose Ud = w c — Xa and consider the Hamiltonian 



H(t) = X a i^j flta + £ o(t)(a + a f ). 



(24) 



Let us denote the computational basis states of N qubits by \ii, . . . , in), with <r 7 j\ij) = ij\ij) and ij £ { — 1, 1}. Under 

the action of the Hamiltonian (24), the initial product state |O)J'0) Ar with \ip) N = i N c h *jv • • • > *iv) evolves 

into 



Jo 



{T yX(r-t)A n iN (25) 



where A lu ... M = 

For simplicity, we consider a square pulse with envelope £o(i) = (qO(T — t). The amplitude of the coherent state 
associated with the qubit component . . . ,In) at time T becomes 

a T = -ie Te 5 smc (XA il ,..., ilt T/2) , (26) 

where sinc(cc) = sin(x)/x. Thus there is both a qubit-state dependent rotation and a suppression of the amplitude of 
the coherent state. Since — N < A^^„ t i N < N, if NXT <C 1, we can expand and obtain 

a T = -ie Q T - e T 2 X ^^ + 0{{NxT) 2 ). (27) 

Thus the first-order correction to the displacement of the cavity state is a qubit-state dependent rotation. The change 
in amplitude is second order in NXT . This has important consequences for the encoding protocol, since it allows for 
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a partial compensation of the finite pulse-length effect of the displacement simply by shortening the duration of the 
subsequent free evolution part. Retaining only the leading order term, the state in (25) factorizes and we have 



\r/>{T)) « \-ie Q T)\i>) N . 



(28) 



which corresponds to a cavity displaced by a = —ie^T, while the qubits were left unchanged. In terms of the drive 
strength and final coherent state amplitude, this is valid when \a\ -C eo/(NX). Because — N < A il ^,^ iN < N, the 
maximum amplitude difference is given to first order by eoT 2 NX = \a\NXT. We may then define a fidelity of the 
displacement operation as the absolute value of the overlap between an ideally displaced coherent state and a state 
with maximal phase difference. To leading order, this is found to be 



cxp 



\{\a\NTxf 



(29) 



Thus we see that to achieve high fidelity, the average number of photons in the cavity n = \a\ 2 should satisfy 
y/fiNXT <C 8. Fig. 4 shows the numerically computed fidelity to the target displaced state obtained from the initial 
state \0) c \gggg) as a function of the square pulse duration T. For this simulation we took X/(2ir) = 10 MHz, a = 2.5. 




0.2 0.4 0.6 0.8 1.0 



N X T/tt 

FIG. 4. Fidelity of un-conditional displacement for the initial state \0) Jgggg) N=4 as a function of square pulse duration T. The 
dashed (red) curve shows F = exp[— (|a|A r xT) 2 /8]. We took a dispersive shift x/(27r) = 10 MHz, and a cavity displacement 
a = 2.5. 



Unconditional qubit selective n-pulses 



In order to perform subset-selective parity measurements, one needs to perform unconditional 7r-pulscs on subsets 
of qubits as explained in the main text. The qubits can be driven simultaneously by using a polychromatic drive. 
Following the discussion in the main text we choose the drive frequencies to be ujf = + 2hX and thus consider the 
following Hamiltonian in the rotating frame 

H(t) = X (j2 at) ota + £ ^ (; """ X V + h.c.) . (30) 
Vi=i / jesi, 

Note that here we have neglected the off-resonant driving of the qubits in Su- In order for this to be valid, we must 
require the frequency detuning between any two qubits to be large compared with 2\ffiX. Moving to a frame rotating 
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with the qubits in S M at frequency nX, we can write the Hamiltonian as (we use the same symbol H to denote the 
Hamiltonian in the new frame) 

H = X \ E ^Wa + xf ^)(a f a-n)+ £ . (31) 

Clearly the first term on the rhs commutes with the remaining terms and leads to a "free" evolution of the qubits in 
Sm and of the cavity given by Us M (t) = exp(— itXa^a ^2 ieSM erf)- The action of the last two terms on each qubit in 
is described by a single qubit Hamiltonian of the form 

H 1 =X(T z (a^a-n) + ^-(T x . (32) 

Defining (3(t) = J* ^rp-dr, the associated evolution operator is 

C/l(t) = e -*/5( t )°- X e -^( ata - fi )/o rfr{ CT *-[co S 2 (/3(T))- S in 2 (/J(r))]+2 -!' S in( / 9(r))co S (/3(r))}^ (33) 

For X = this reduces to a simple cr x rotation, the angle of which is determined by the integral of the pulse envelope. 
The total evolution operator, including the contribution from the first term on the rhs of (31), is then 

U(t) = U S {t)e~ imT ' 3eS M ^ e -^( at °-™)E )e sS ; / t dr{ CT j[cos 2 (/3(r))- S in 2 (/3(r))]+2 CT J S in(/3(r))cos(/J(r))} ^ 

For a square pulse envelope = Vl^0{T — t) with T = 7r/Oo such that (3(T) = ir/2 it is easy to show that 

I p 2 ^3 £S M 3 P n V ' A "^ S M 3 



f/(T) =C/ Sm [jf) e ~^ 3es " 3 e n ^"^z., 6S& ",_ (35) 

We are interested in the action of this operator on a coherent state or superposition of coherent states, for which the 
photon number is Poisson distributed around n with variance equal to n. Hence if 2^/fiX <C we have 



U(T)nU SM (T) 



H) M II 



/2y/ExT\ 



\ TT J 

If furthermore, we require that 2XMir < f2 , then J7 Sjf « 1 + 0(2XMT) so that 



t/(T) » H) M J] trf + O (^^,2xMt\ 



(36) 



(37) 



Fig. 5 shows the fidelity as a function of the square pulse duration T, of an unconditional two-qubit tt rotation 
acting on an initial state |a) c |^>} 4 with \i))) A = y/l/'S(\eegg) + \eggg) — \egge)), a = 2.5 and for X/(2tt) = 10 MHz. 

STABILIZER PUMPING 

In the main text we show how to measure a stabilizer by first encoding its eigenvalues onto the cavity state and 
then reading out the latter either directly or via the ancilla qubit and the readout cavity. Alternatively one can also 
implement conditional jump operators of the form ar'fPg M where -Pfj^ = (1 ± Zs M )/2 with Zs M = Y[ies M °"i an< ^ 
i £ Sm . Such operators can be used to digitally simulate an artificial bath that pumps the system onto the stabilizer 
subspace as first proposed for trapped ions in [14] . The dissipative evolution under the action of such jump operators 
can be implemented with the protocol shown in Fig. 6. As for a parity measurement, one first encodes the eigenvalues 
of the parity operator onto the cavity state and applies and unconditional displacement to obtain the state (see Eq. (5) 
of the main text) 

| tt) = 1 2& N )P+ M |V) W + |0) P Sm \4>) n . (38) 

Instead of mapping the cavity state onto an ancilla as in the measurement protocol, one then applies a rotation by 
an angle 9 around x of qubit i € Sm conditioned on the cavity being in the vacuum state. Note that this does not 
require an ancilla qubit. Because {erf, Z$ M }+ = for i € Sm, the state becomes 

P+ M (|^) JV |2a) c + i S in(e/2K|^) iV |0) c ) + co S ^/2)P SM |^|0> c . (39) 
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nxT/n 

FIG. 5. Fidelity of un-conditional tt rotation of the first two qubits of the initial state |o) c |V , )4 as a function of square pulse 
duration T. The cavity displacement is a = 2.5 and the dispersive shift is x/(27r) = 10 MHz. 

By resetting the cavity to the vacuum (via an ancilla qubit or a cavity Q-switch [26]) and repeating the procedure 
one generates the dissipative (non-unitary) evolution [14] 



® Pn 



E o (0)p N El(0) + E^(0)p N E^(0) 



(40) 



t- cos(0/2)Pg M and E { x t] (0) = isin(0/2)of P^. Thus with probability sin^(6>/2), a -1 eigen- 



with Eq{6) = PcJ 

state of Zs M is converted into a +1 eigenstate. For 6 — tt, this "stabilizer pumping" occurs with unit probability. 
Generalization to arbitrary weight-M Pauli-operators can be achieved with suitable single qubit rotations. 



|0) c -\D^ 
1 



N 



it 



D a 



n. 



reset|0) c 



Y° 



FIG. 6. Quantum circuit diagram for stabilizer pumping. D a denotes an unconditional displacement by a. Ti/ 2 represents a 
free evolution of duration Ti/ 2 = t/4x- H-Sf, represents a (parallel) multi-qubit 7r-pulse on qubits in S%[. Xf e s M ,e represents 
a conditional rotation of qubit i by 6 around the a;-axis, conditioned on the cavity being in the vacuum, reset |0) c represents 
the (dissipative) cavity reset operation. 



DECOHERENCE 

Decoherence due to photon losses and qubit decay and dephasing are treated within the Lindblad master equation 
formalism [33]. The master equation we solve numerically is 



[H, p] + KD[a]p + ( -fVWfiP + r-2?[o-r]p + +T+D[tr+], 



(41) 
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where T>[L]p = (2L/>L^ — L'Lp — pL J L)/2. Here k is the single photon decay rate of the cavity, T v is the phase 
relaxation (pure dephasing) rate of the qubits and the qubit relaxation rate is given by 1/T\ = T_ We consider 

the zero-temperature limit T + = 0. 



LOGICAL QUBIT STATE TOMOGRAPHY 



In Fig. 3 of the main text, we show the state fidelity of the encoded logical qubit state to the target state 
exp[— z(7r/8)X]|0) as a function of the dispersive shift X and cavity nonlincarity K. In the left panel of Fig. 7, 
we compare the resulting four-qubit state with the ideal target state in the Pauli-bar representation of the reduced 
density matrix (tracing out the cavity), for K/(2n) = 80 KHz, X/(2tt) = 5 MHz and k/(2tt) = 10 KHz. The fidelity 
to the target state is 94%. The right panel of Fig. 7 shows the Bloch-sphere representation of the prepared logical 
qubit state. From the latter, we sec that the infidelity is due to both a loss of purity and a deviation from the ideal 
rotation. 
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FIG. 7. (Color online) Left panel: Pauli-bar representation of the prepared four qubit state (full (blue) bars) and of the target 
state exp[-i(7r/8)X]|0) (lighter (red) bars). Here P 4 = {1, cr x , cr v , a z }® 4 is the Pauli- group of order four. Right panel: 
Bloch-sphere representation of the prepared logical qubit state (full (blue) vector) and corresponding target state (full (red) 
vector) . 



DISCUSSION & OUTLOOK 

We conclude these notes by an informal discussion of a few possible applications of our protocol to illustrate its 
versatility. 



Implementation of the toric code 

Because we can measure Pauli operators on arbitrary subsets of qubits, it becomes possible to implement "in 
software" non-trivial topologies such as required e.g. by the toric code [12], in which 2g logical qubits can be encoded 
into the ground state stabilizer subspace of a qubit manifold with genus q. Let us consider the case q — 1 of a ring 
torus as depicted in Fig. 8. The stabilizers of the 27V-qubit toric code for example are the N plaquettc Pi and N star 
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FIG. 8. Left panel. 8-qubit toric code. The qubits are represented by read beads. Right panel. Corresponding planar 
representation. The latter is easier to implement experimentally and the toric boundary conditions are imposed when defining 
the stabilizers. C t indicate the incontractible loops used to define the logical qubit operators (see text). 



operators Si defined as 

Pi= n s t = n ( 42 ) 

J tzplaquctte^ j€star^ 

Crucially, because of the toric boundary conditions it holds true that riigiattice ^ = ELeiattice Si = 1. These two 
constraints then lead to a four-fold degeneracy (2 2fl -fold degeneracy in general) of the ground state manifold. A 
ground state of the toric code can be prepared by starting from the fully polarized state \g, . . . , g) of the 2N qubits in 
the Z basis and sequentially measuring all the star operators and applying a corrective single qubit Z operation on a 
qubit in any star with eigenvalue —1. Fig. 8 illustrates a minimal 8-qubit toric code and its planar implementation. 
The stabilizers of this code are explicitly: 

Pi = Z\Z±Z-jZ% Si = X1X2X5XJ 

P2 = Z^Z^Z'jZ^ S2 = X^X/iX^Xy 

P 3 = ZiZ4Z 5 Z e S3 = XiX 2 X e Xg, 

P3 = Z2Z3Z5Z6 S3 = X^X^XqXs (43) 

Here Zi = erf and X; = erf. Our protocol can also be used to measure the logical Pauli operators of the toric code, 
which arc defined in terms of the incontractible X or Z loop operators winding the torus in either vertical or horizontal 
directions (see Fig. 8). Explicitly 

z (i) = n < z m - n « 

X {1) = J] of X i2) = II at- (44) 



Preparation of cluster states &: measurement based quantum computation 

Cluster states are the central resource of measurement based quantum computation [11]. A cluster state is a 
stabilizer state defined on a graph G, where each vertex i is associated with a stabilizer generator of the form 
Si = Xi rijenn(i) wnere nn(i) denotes nearest neighbor vertices to vertex i. By applying a single qubit Hadamard 
rotation on qubit i at the beginning and at the end, we can use our encoding procedure to measure Si and hence prepare 
arbitrary cluster states. To perform a quantum computation one requires in addition the ability to perform singlc-qubit 
measurements. Because single qubit operators are weight-one Pauli operators, we can again use our measurement 
protocol to perform such measurements and hence perform one-way measurement based quantum computation. In 
practice these single-qubit measurements can be done either as explained in the main text, by applying parallel echo 
7r-pulses on all the qubits except the one being measured, or alternative by first performing a parity measurement of 
all the (appropriately rotated) qubits (i.e. without any echo 7r-pulses) yielding P t ot G {il} an d then repeating the 
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measurement, echoing-out only the qubit that is to be measured yielding P C omp € {±1}- The parity of the qubit to 
be measured is then given by Pi = Ptot/-Fcomp- The advantage of the latter approach is a drastic reduction in the 
number of necessary 7r-pulses (from N — 1 to 1 where N is the total number of qubits). This reduction applies in 
general for subset parity measurements when the total number of qubits N is larger than twice the number of qubits 
in the subset M. 

Conclusion 

The proposed protocols offer a direct route to implement stabilizer quantum error correction and to "digitally 
simulate " models of topological quantum computation with superconducting qubits in the ultra-strong dispersive 
regime. From an experimental point of view, the main challenge presumably lies in engineering a system of multiple 
qubits with long coherence times and with sufficient separation of time scales to validate the approximations made in 
our derivations. A proof of principle experiment seems within reach of current technology. 

As an outlook for future research, an interesting question is whether our approach can be extended to measure all 
stabilizers of a given code in parallel rather than sequentially. This should be possible in principle as these operators 
all commute with each other. A naive extension to encode an entire syndrome (collection of N stabilizers, each of 
which can have values ±1) into 2 Ar -component cavity cat-states, would seem to require an unfortunate exponential 
increase in the required cat-size. 

I 



